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Abstract 

We show that the failure of the real-space RG method in the ID tight-binding model is not 
intrinsic to the method as considered so far but depends on the choice of boundary conditions. For 
fixed BC's the failure does happen. For free BC's we present a new analytical block RG-method 
which gives the exact ground state of the model and the correct 1/iV^-law for the energy of the first 
excited state in the large iV(size)-limit. We also give a reconstruction method for the wave-functions 
of the excited states. 
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1 Introduction 



Real space Renormalization Group (KG) methods originated from the study of the Kondo problem 
by Wilson Q. It was clear from the beginning that one could not hope to achieve the accuracy 
Wilson obtained for the Kondo problem when dealing with more complicated many-body quantum 
Hamiltonians. The key difference is that in the Kondo model there exists a recursion relation for 
Hamiltonians at each step of the RG-elimination of degrees of freedom. Squematically, 

Hn+1 = Hi\j + hopping boundary term (1) 



Hj,+i = R{Hn) (2) 

where R is the RG transformation. 

The existence of such recursion relation facilitates enormously the work, but as it happens it is specific 
of impurity problems. 

From the numerical point of view, the Block Renormalization Group (BRG) procedure proved to 
be not fully reliable in the past particularly in comparison with other numerical approaches, such as 
the Quantum MonteCarlo method which were being developed at the same time. This was one of the 
reasons why the BRG methods remained undeveloped during the '80's until the begining of the '90's 
when they are making a comeback as one of the most powerful numerical tools when dealing with zero 
temperature properties of many-body systems, a situation where the Quantum MonteCarlo methods 
happen to be particularly badly behaved as far as fermionic systems is concerned 

As it happens, the BRG gives a good qualitative picture of many properties exhibited by quantum 
lattice Hamiltonians: Fixed points, RG-flow, phases of the system etc. as well as good quantitative 
results for some properties such as ground state energy and others Q, For a review on the 

Block RG method see |^ and chapter 11 of reference |^]. However in some important instances the 
BRG method is off the correct values of critical exponents by a sensible amount. 

The first advance in trying to understand the sometimes bad numerical performance of the BRG 
methods came in the understanding of the effect of boundary conditions (BC) on the standard RG 
procedure 

White and Noack pointed out that the standard BRG approach of neglecting all connections 
to the neighbouring blocks during the diagonalization of the block Hamiltonian Hb introduces large 
errors which cannot be corrected by any reasonable increase in the number of states kept. Moreover, 
in order to isolate the origin of this problem they study an extremely simple model: a free particle in 
a ID lattice. As a matter of fact, it was Wilson ||l^ who pointed out the importance of understanding 
real-space RG in the context of this simple tight-binding model where the standard BRG clearly fails 
as we are going to show. 

The reason for this failure can be traced back to the importance of the boundary conditions in diagonal- 
izing the states of a given block Hamiltonian Hb in which the lattice is decomposed into. Notice that 
in this fashion we are isolating a given block from the rest of the lattice and this applies a particular 
BC to the block. However, the block is not truly isolated! A statement which is the more relevant 
the more strongly correlated is the system under consideration. Thus, if the rest of the lattice were 
there it would apply different BC's to the boundaries of the block. This in turn makes the standard 
block-diagonalization conceptually not faithfully suited to account for the interaction with the rest of 
the lattice. 

Once the origin of the problem is brought about the solution is also apparent: devise a method 
to change the boundary conditions in the block in order to mimick the interaction with the rest of 
the lattice. This is called the Combination of Boundary Conditions (CBC) method which yields very 
good numerical results. This method has not yet been generalized to interacting systems. However in 
reference [11 1 an alternative approach is proposed under the name of Density Matrix Renormalization 
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Group (DMRG) which apphes to more general situations and also produces quite accurate results. 
Despite of the impresive numerical accuracy of these methods there is not a clear understanding of why 
they perform so well or whether they would be valid in dimensions higher than one. 



2 Boundary Conditions and Real-Space RG 

In this letter we want to reconsider again the role of BC's in the real space RG method for the case of 
a single-particle problem in a box. The continuum version of this Hamiltonian is simply H = 
We shall consider open chains with two types of BC's at the ends: 

Fixed BC's: ^(0) = i^{L) = 



Free BC's: ^(0) = |^(L) = 
ox ox 

The lattice version of H for each type of BC's is given as follows: 



(3) 
(4) 
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The only difference between Hpixed and Hpree appear at the first and last diagonal entry (2 ^ 1). The 
exact solution of (^) and is very well-known and we give it for completeness: 

Fixed BC's: 



N + 1 



4sm^(— ^ - 



1,2, 



n = 0,1, 



Free BC's: Vn(i) = <^ cos ^(j - i), = Asm\^) 



(7) 
(8) 



.A^-1. 



where the N^^s are normalization constants and is the number of sites of the chain. 

Before getting into the problem of the renormalization of these Hamiltonians, it is worth to pointing 
out another physical realization of Hpree- Consider the following Hamiltonian for the anisotropic 
Heisenberg model in an open chain of N sites, 



N-l 



Ha = -IJ: (ajaj+i + + A(a|a|+i - 1)) 



(9) 



where A = -|-1 (A = —1) corresponds to the Ferromagnetic case (Antiferromagnetic case). 
When A > 1, the ground state of is the Ferromagnetic state \F) given by: 
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|F) = 1 TT...T) 



(10) 



Ha\F) = 



:ii) 



Now, let us denote by \x) the state which is obtained from the Ferromagnetic state \F) when one spin 
is overturned at the position x of the chain. The subspace spanned by the \x) states is the 1 magnon 
subspace of the Ferromagnetic model. It is readly checked that the restriction of the Ferromagnetic 
Heisenberg Hamiltonian takes the following form: 



HMagnon{^) 



( 1 -1 

-1 2A -1 
-1 2A 



V 



2A -1 
-1 1 / 



(12) 



If A = +1 it precisely coincides with Hpree- 



HMagnoni^ — +1) — Hpree (13) 

Thus we arrive at the following mapping: A simple magnon above a ferromagnetic background satisfies 
Free EC's. We shall come back to this connection at the conclusions for it opens the door to the 
introduction of interactions in the model. 

Now let us get to the problem of renormalizing the tight-binding Hamiltonians (^)-(|6|). 

In figures 1 and 2 we show the ground state and first excited states of the chain with fixed and free 
EC's. 

It is clear from Fig.l that a standard Block RG method is not appropiate to study the ground state 
of fixed BC's since this state is non-homogeneous while the block truncation does not take into account 
this fact. 

Each piece of the ground state within each block satisfies BC's which vary from block to block. This is 
the motivation of reference to consider different BC's in the block method, yielding quite accurate 
results as can be seen from Table 1 for comparison of the CBC method and the standard method. 
We observe that the standard RG method performs rather poorly as compared to the CBC method 
which yields quite the exact results. 

The other alternative to the CBC method is the Density Matrix RG method which can be phrased by 
saying that the rest of the chain produces on every block the appropiate BC's to be applied to its ends, 
and it has the virtue that can be generalized to other models, something which is not the case as for 
the CBC method. 

On the other hand, the ground state of Hpree is an homogeneous state (see Fig. 2) which in turn suggests 
that a standard RG analysis may work for this type of BC's. We shall show that this is indeed the 
case if the RG procedure is properly defined. The key of our RG-prescription is to notice that Hp^^gg 
has a geometrical meaning: Hpree is the incidence matrix of the graph in Fig. 3, and it is called minus 
the discrete laplacian —A of that graph. Notice that Hpi^f.^ has not such geometrical interpretation, 
in fact, it concides with the Dynkin diagram of the algebra A^. 

Based on this observation the Kadanoff blocking is nothing but the breaking of the graph into N/ug 
disconnected graphs of rig sites each. We shall choose = 3 in our later computation as shown in Fig. 
4. 

The previous geometrical interpretation of Hpree suggests that we choose the block Hamiltonian 
Hb to be the incidence matrix of the disconnected graph in Fig. 4, namely, 
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(14) 



in which case the interblock Hamiltonian Hqb which describes the interaction between blocks becomes: 



/ 



Hbb 



\ 



-1 
1 



(15) 



V ■•• / 

Hbb in turn also coincides with the incidence matrix of a graph which contains the links missing in 
Fig.4 (RHS) which connects consecutive blocks. In a few words: our RG-prescription introduces free 
BC's at the ends of every block. This condition fixes uniquely the breaking of Hp^ee into the sum 
Hb + Hbb- This is the choice we make. It should be emphasized that the splitting of Hpree into two 
parts Hb + Hbb is by no means unique, so that different choices may lead to very different results. 
Prior to any computation we notice that the previous RG-prescription should lead to an exact value of 
the ground state energy, for the ground state of each block is again a constant function. The question 
is therefore to what extent our method is capable of describing the excited states. We shall concentrate 
ourselves to the first excited state since computations can be carried OTit analytically. 

First of all we diagonalize Hb within each block of 3 sites, keeping only the ground state V'o^^ a^nd 
the first excited state if)^^ (3 — 2 truncation) . The superscript denotes the initial step in the truncation 
method. In the standard RG method we would choose and V'j^^ as the orthonormal basis for the 
truncated Hilbert space and obtain the effective Hamiltonian H'^ and H'^^. In our case it is convenient 
to express these effective Hamiltonians in a basis expanded by the following linear combination: 



1 
72 



(0) 



(16) 
(17) 



which are also an orthonormal basis of the truncated Hilbert space. 
In this basis the truncation of Hb reads as follows, 



Hb — *■ HL 



( ^ 








\ 
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A 
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V 






I 



(18) 
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with e taking on the value e'-'^^ = 1 in the initial step of the RG-method, which is the energy of the 
state V'l''''- 

The truncation of Hbb is more complicated, the result being: 



Hbb — > H'r 



BB 



( B C 
D + B C 
C* D + B 



D + B C 
C* D 



(19) 
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\ and 6(°) 



+ i in the initial step of the RG-method. 



with a 6 taking on the values a^^^ = 

The nice feature about the basis ([l^-(p!7|) is that all rows and columns of (18) and (|l^) add up 
to zero, just like the original Hamiltonians ( [T4|) and ([Tsl), implying that the constant vector is an 
eigenvector with zero eigenvalue of the renormalized Hamiltonian! 

We shall call H]y/^{e,a,h) the sum of the Hamiltonians (|T8| ) and (|l9|) for generic values of e, a and 
6. Next step in our RG-procedure is to form blocks of 4 states of the new Hamiltonian Hj\f^-^{e,a,h) 

and truncating to the two lowest -00^^ and tp^^ energy states within each 4-block (4^2 truncation). 
The reason for this change in the number of sites per block (from 3 to 4) is motivated by the form of 
H'^Q in (|l9|) and the fact that if we try to make a second step in the RG-method with 3-blocks the 
method is doomed to failure because the constant state of Fig. 2 would no longer be the ground state. 
Fortunately enough, with 4-blocks if we define new states "0+ ^ and tp^^ in the same form as we did in 
Eq.(|5|), we obtain that the new effective Hamiltonian is obtained by a redefinition of the parameters, 
namely, 



Hn/^{(-, «, ^) — ' H,y/Q{e', a', h') 



(20) 



1 + a2 ^ 52 _ ^ 



(21) 



1 



6' 



2^/2 
1 

2^/2 



a + b- 



a + h + 



a{a^ - 362 + a) + 



v/A(A + a2 - 62) 



a{a' - 362 + A) + I 
v/A(A + a2^62) 



(22) 



(23) 



A = W(a2-62)2 + (l-2a6)2 



(24) 



In this fashion, the constant state of Fig. 2 is again the ground state of the model and moreover, 
upon iteration of Eqs.(|20|)-(|2^) there are no level crossing among the excited states. Otherwise stated 
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this means that the level structure of the block Hamiltonian Hb is preserved under the action of our 
BRG-method based upon the reduction from 4 to 2 states. 

The energy Ei{N) of the first excited state of a chain with = 3 x 2™ sites can be obtained 
iterating m times Eqs.pO|)-(p^): 

Si(iV = 3x 2™) = e(") (25) 

The initial data are given by: 

e(o) = i, «(0)= 1 _1 ^(0)^ 1 +1 (26) 
^ 2' ^ 2 ^ ^ 

In Table 2 we give our results for small and large number of iterations. For low values of the 
deviation of e^"^^ with respect to the exact result is small. Recall that we are only keeping two states in 
our RG-procedure, and that the ground state energy is exactly zero by construction!. But what is more 
interesting about these results in Table 2 is that we are able to obtain the correct size dependence, i.e., 
of e^"*). As a matter of fact, the energy of the first excited state behaves for large A^ as (^): 

E[^---'\N) ^ Cec.act/N^ With Ce^act = (27) 

while our BRG-method gives, 

e[^R(^\n) cbrg/N^ with CBRG = 12.6751 (28) 

The achievement of the 1/A^^-law is a remarkable result which in turn allows us to match the correct 
order of magnitude of the energy. For instance, for 10 iterations our RG-method with 2 states kept 
gives the energy of the order of 10~^, which is precisely the same order of magnitude as for the CBC 
method (see Table 1) but with 8 states kept in the case of Fixed BC's. Recall that the standard BRG 
performs as bad as a 10~^ order of magnitude. 



3 Wave- Function Reconstruction 

We may also wonder whether we are able to make a reasonable picture of the first excited state wave- 
function based upon our BRG-procedure when compared with the exact form depicted in Fig. 2. As we 
are working with a real-space realization of the renormalization group method, this is something we have 
at hand. To do this we need to perform a "reconstruction" of the wave-function. This reconstruction 
amounts to plot the form of our aproximate wave-function in each and every of the 3-blocks out of 
the 2""+^ in which the original chain is decomposed into under the BRG-procedure. Recall that in the 
initial step we started out with blocks of 3 states keeping the two lowest states Vo'^^ and (3 — > 2 
truncation). In the next step we make blocks of 4 states keeping the two lowest states V'o^^ and tp^'^ 
(4 — > 2 truncation) and then we perfom the iteration procedure over and over. As a result of this 
procedure we may express the two lowest wave functions of the m + 1-th step in terms of those of the 
previous m-th step by means of the following matricial form: 



; (m+l) 




-Ctr. 




(29) 

R 

where the LHS of Eq. ( p9| ) represents the wave function of 3 x 2"*^"^ sites while in the RHS we have a 
left-wave-function of 3 x 2™ sites and another right-wave-function of 3 x 2™ sites, so that everything 
squares. The parameters appearing in Eq. (|29| ) turn out to be given by: 
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- 2amhm) + (q^ - fe^ + a, 
2v/A„(A™ + a^-6^) 

- 2a. J)^) - (o^ - + A^ 
2v/A„(A^ + a^-6^) 



(30) 
(31) 



with Am as in Eq.(|2^ 
We may recast Eq. 



where 



Their initial values are ao = 1/vlO and /?o = 3/vlO. 
in more compact form by writing: 



(32) 



Ijr, 
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' 1 
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Til 




Pm 



(33) 



(34) 



We may call Eq.(|3^) the reconstruction equation. This is the master equation that when iterated 
"downwards" (reconstruction) allows us to obtain the picture of our approximate BRG-wave-function 
corresponding to every and each block of 3 sites of the 2*"+^ blocks in which the chain is decomposed 
into. At the end of the iteration procedure we end up with expressions for the values of the 3-sites 
wave- functions in terms of the initial two lowest states Tp^^^^ and tl^f'^. The first one is a constant 
function while the second is a straight line of negative slope. Thus, these two states turn out to be 
the building blocks of our BRG-procedure. As a matter of illustration, we present in Fig. 5 the plot for 
the reconstruction of the first excited wave-function for a chain of = 48 (m=3) sites long. We can 
observe that our BRG-wave-function has the shape of a broken line around the exact form of the wave 
function. We have checked that this broken shape remains when increasing the length of the chain so 
that its nature must be due to the fact that we are only keeping 2 lowest states in our procedure. It 
is expectable that upon increasing this number the shape of the aproximate wave function must be 
smoothed out. Nevertheless, the important point is that our BRG-method preserves the number of 
nodes of the wave-function which one of the defining features of the exact result. 

When using the reconstruction equation to obtain the wave function we may use a binary code based 
upon the labels L (left) and R (right) to keep track of the different 3-sites blocks which make up the 
chain. Thus, in one dimension the RG-blocks are in a one-to-one correspondence with a binary numerical 
system. In general, for other dimensions we may state squematically the following correspondence: 



BRG-prescription 



"Number System" 



4 Conclusions and Outlook 

In this letter we have shown that it is possible to perform a successful real-space block RG treatment of 
the ID tight-binding model. In doing so we have elucidated the role played by the boundary conditions 
in the real-space RG method for we have shown that the failure of the BRG-method pointed out in 
1^] can be traced back to a particular choice of boundary conditions: if we choose free BC's instead of 
fixed BC's the ground state turns out to be homogeneous, a property which is necessary in order to 
have a suitable implementation of a blocking procedure. Otherwise, for fixed BC's the blocks of the 
RG-method cannot have fixed BC's as well for in that case one introduces anomalous behaviour in the 
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ground state wave-funcion which ruins the method. Thus, the failure of the BRG-method accepted so 
far is not intrinsic to the method but depends on the choice of the boundary conditions. 
We have devised a BRG-method which always yields the correct ground state of the Hamiltonian with 
free EC's for open chains of any length. As for the excited states, our method is able to reproduce 
the correct 1/A^^-law exhibited by the exact first excited state (|27|)-(|28|) in the large A^- limit with the 
correct order of magnitud. Since we only keep 2 states - to make computations analytical - at every 
step of our truncation procedure, this means that we have been able to capture the essential physics of 
the model in the most economic way. For this reason we believe that our method should also be able 
to account for the correct properties of the rest of excited states of the model. 

We have given another physical interpretation for the tight-binding Hamiltonian with free EC's as the 
restriction of the Ferromagnetic Heisenberg model to the subspace of one-magnon solutions (IT^ ) . In the 
light of this mapping we may interpret the 1/A^^-law of our ERG-method as in good agreement with the 
correct quadratic dispersion relation exhibited by the Ferromagnetic magnons E{k) ^ k"^ ^ l/N"^ for 
low wave number k. Thus, from this point of view, our RG-procedure is aplicable to the Ferromagnetic 
regime of the Heisenberg model and we may wonder whether we may extend it to the Antiferromagnetic 
regime and obtain the linear dispersion relation E{k) ~ A; ~ 1/A^. 

Our ERG-method also gives a good qualitative picture of the ground state and first excited states 
wave- functions in coordinate space as shown in Fig. 5. In fact, the number of nodes is preserved by the 
ERG-procedure. This good performance is due to a reasonable choice of the 2 lowest energy states kept 
in the initial stage of the method. The reconstruction of the wave-function is based upon the master 
equation (^2|) which can be readly generalized to account for a bigger number of states kept. 
Finally, we may address yet another generalization of the present work by considering a Hamiltonian 
HpTeeil) as in Eq.(^) but depending upon a complex parameter q. This parameter enters in the 
matrix expression of the Hamiltonian by replacing the diagonal (1, 2, 2, . . . , 2, 1) of Eq.(^ by (g, q + 
q~^,q + q~^, . . . ,q + q^^,q^^) while the off-diagonal terms remain as -1. This corresponds to a one- 
parameter deformation of the boundary conditions in which the free case is recovered as the q ^ 1 
limit. Furthermore, it is more interesting to consider this extension from the point of view of the 
Heisenberg-model connection for in this case the Hamiltonian comes with a boundary term derpending 
on the q parameter: 

The important point about this Hamiltonian is that it is invariant under what is called the quantum 
group SUq{2). We have shown in that this symmetry plays an important role in constructing an 
improved version of the standard ERG-procedure. We call this new method a g-ERG-method and for 
instance it is able to obtain correctly the line of critical XXZ models, unlike the standard ERG-method. 
This quantum symmetry must show up in the free case when performing the previous g-deformation. 



e-mail addresses: mardel@fis.ucm.es (M.A. M.-D.); sierra@cc.csic.es (G.S.) 
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Table captions 

Table 1 : Exact, Standard RG and CBC Values of Low Lying States for the ID Tight-Binding Model 
for a chain of iV = 2048 sites. 

Table 2 : Exact and new standard RG values of the first excited state for the ID Tight-Binding 
Model. 
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Figure captions 

Figure 1: Ground state ipo and first excited state ipi for the Hamiltonian Hpixed with fixed BC's. 
Figure 2 : Ground state V'o and first excited state for the Hamiltonian Hpree with fixed BC's. 
Figure 3 : The onc-dimensional open chain as a graph. 
Figures 4 : Lattice decomposition into blocks of Ug sites (n^ = 3). 

Figure 5 : The reconstruction of the first excited wave-function for a chain of A?" = 48 (m=3) sites 
long. 
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Energies 


Exact 


Standard RG 


CBC 


Eo 


2.351 X 10"*^ 


1.9207 X 10-^ 


2.351 X 10-** 


El 


9.403 X lO"*^ 


1.9209 X 10"^ 


9.403 X 10-*^ 


E2 


2.116 X 10-^ 


1.9214 X 10-^ 


2.116 X 10-^ 



Table 1: Exact, Standard RG and CBC Values of Low Lying States for the ID Tight-Binding Model 
for a chain of N = 2048 sites. 
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12.6751/Af^ 



Table 2: Exact and new standard RG values of the first excited state for the ID Tight-Binding Model. 
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